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The experimental determination of the properties of the newly discovered boson at the Large 
Hadron Collider is currently the most crucial task in high energy physics. We show how information 
about the spin, parity, and, more generally, the tensor structure of the boson couplings can be 
obtained by studying angular and mass distributions of events in which the resonance decays to 
pairs of gauge bosons, ZZ, WW, and 77. A complete Monte Carlo simulation of the process pp — > 
, X — > VV 4/ is performed and verified by comparing it to an analytic calculation of the decay 

amplitudes X — > VV — > 4/. Our studies account for all spin correlations and include general 
couplings of a spin J = 0, 1,2 resonance to Standard Model particles. We also discuss how to use 
angular and mass distributions of the resonance decay products for optimal background rejection. 
It is shown that by the end of the 8 TeV run of the LHC, it might be possible to separate extreme 
hypotheses of the spin and parity of the new boson with a confidence level of 99% or better for a 
wide range of models. We briefly discuss the feasibility of testing scenarios where the resonances is 
not a parity eigenstate. 
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> " I. INTRODUCTION 

00' 

The discovery of the new boson [l], Q at the LHC, which is further corroborated by the strong evidence from the 
Tevatron is the culmination of the hunt for the elusive Higgs boson. Three primary decay channels 1 X — s- ZZ, 
[ WW, and 77 were observed experimentally by the CMS and ATLAS collaborations. However, not much is currently 
known about detailed properties of the new boson beyond its mass, nix ~ 125 GeV, although some information 
\ can be reasoned from data. We know that the width of the new particle is consistent with being smaller than the 
experimental resolution of about a GeV. We also know that, as a consequence of the Landau- Yang theorem d, [f|, 
the new boson cannot have spin one because it decays to two on-shcll photons. Finally, we know that the relative 
decay branching fractions and production cross-sections of the new particle are generally consistent with the Standard 
Model (SM) Higgs boson hypothesis although current accuracy of experimental measurements does not allow for 
an unambiguous conclusion. 

Since the new boson interacts with massive gauge bosons, we expect it to play some role in electroweak symmetry 
breaking. However, this needs to be verified by direct measurements of its properties. In particular, it is important to 
experimentally study the tensor structure of couplings of the new boson to SM fields and its SU{2) x £7(1) quantum 
numbers (if any), avoiding theoretical prejudice. For example, we may wonder if the relatively strong interaction of 
the new particle with electroweak gauge bosons already observed implies that this new boson is not a pseudoscalar. 
One may argue that this is the case because a pseudoscalar must interact with gauge bosons by means of higher- 
dimensional operator whose significant contributions to X — > VV would imply low scale physics beyond the SM which 
should have already been observed experimentally. Since no beyond the SM physics has been observed at the LHC, 
the scale of new physics cannot be low and it is tempting to conclude that the pseudoscalar nature of the new boson is 
excluded 2 . While such arguments are appealing and may in fact be valid, it is important to test them experimentally, 
especially when such tests are within reach. In fact, as we show in this paper, it is entirely possible to achieve that 



1 Throughout this paper, we will use a uniform notation for both on-shell and off-shell massive gauge bosons. 

2 Even stronger arguments about XVV coupling by means of higher-dimensional operators are possible if one assumes that X is a singlet 
under electroweak SU{2) X U(l) 0. 
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with data from the 8 TeV run of the LHC using di-boson final states. Hence, it is realistic to expect that a clear 
profile of the new boson can be established by purely experimental means in a short period of time. 

The determination of the quantum numbers of a Higgs-like particle was discussed in great detail in the literature, 
see Refs. [8T424|. The strategy that we use in this paper is similar to what has already been discussed in Ref. [20| . 
In that reference we demonstrated that X decaying to two vector bosons provides an excellent channel to study the 
tensor structure of its couplings and outlined the general way to do so. Since Ref. [20( was written before the new 
particle was discovered, its mass and production rates were unknown. As a result, many examples studied in Ref. (20j 
are, by now, of an academic interest. The discovery of the new boson allows us to extend the discussion presented in 
Ref. [20] and arrive at realistic predictions about the prospects for measuring its spin and couplings at the LHC. 

We extend the analysis reported in Ref. [20] in several important ways. First, since the mass of the new resonance 
is m x ~ 125 GeV, at least one of the bosons in the X — > ZZjWW decay is off-shell. Calculations reported in Ref. (2pj 
employed general structures of scattering amplitudes and general angular distributions but did not fully include the 
off-shell kinematics of the vector bosons; we improve on this in the current paper. We also extend those earlier results 
by including WW and 77 final states in the Monte Carlo simulation. As we stressed earlier [2(| , the optimal analysis 
for the new boson discovery and its property measurements requires utilization of the full kinematic information about 
the process. Analysis based on matrix elements or multivariate per-event likelihoods, such as MELA (Matrix Element 
Likelihood Analysis), adopted by CMS [l,[25j], allows for optimal background suppression. The same techniques also 
guarantee the best performance when applied to measurements of the new boson's properties. 

In this paper, we consider the gluon fusion, gg, and quark-antiquark annihilation, qq, production mechanisms. The 
primary production mode of the SM Higgs boson is expected to be gluon fusion. The inclusion of the qq production 
process completes all the possible initial state polarization scenarios for spin-one and spin-two resonance hypotheses 
thus allowing for the most general treatment of kinematics, and inclusion of all relevant spin correlations. We also 
note that weak Vector Boson Fusion (VBF) is expected to account for 7% to the SM Higgs boson production rate. 
Since jet tagging identification would reduce the experimentally observable rate even further, the contribution from 
the VBF topology is at the level of a few percent. As a result, we leave dedicated analysis of the VBF topology as 
well as the analysis of other final states in the decay of the new boson to future work. 

The paper is organized as follows. In Scction|H]we review kinematics in resonance production and decay, expanding 
on our earlier work in Ref. [201 ] and focusing on the case relevant to the observed boson mass mx < 2mz(w)- hi 
Section IIIII we discuss the Monte-Carlo event generator for simulating production and decay of a new boson with 
different hypotheses for spin and tensor structure of interactions, expanding Ref. [20j to include new final states 
covered in this paper. In Section [IV] we discuss the analysis methods. We summarize the results and conclude in 
Section [V] Detailed formulas for angular distributions and some numerical results are given in the Appendix. 



II. KINEMATICS IN RESONANCE PRODUCTION AND DECAY 

Before going into the discussion of how properties of the boson X can be studied, it is interesting to point out that 
the determination of the spin-parity of a resonance through its decays to two gauge bosons, that subsequently decay 
to four leptons, was first attempted more than a half-century ago, with the study of neutral pion decays ir° — > 77 
and 7T° — » 7*7* — » e + e - e + e~. Photon polarization in tt° — > 77 can be used to determine ir° parity |5|, but it is 
more practical to use the orientation of the planes of the Dalitz pairs in the decay n° — > e + e~e + e~ [26]. Further 
developments of these techniques were discussed in Refs. [13, [H[ and additional refinements were suggested in the 
context of Higgs physics in Refs. 8-24] and in the context of -B-physics in Refs. [291432] ]. By analogy, the decay 
X — > ZZ — > At is an excellent channel to measure the spin, parity, and tensor structure of couplings of the new boson 
since the full decay kinematics are experimentally accessible. In the channels X — > WW — > 2l2v and X — > 77 less 
kinematic information is available, but they can complement the measurements of the resonance properties. Other 
final states of ZZ and WW could be considered, but they typically suffer from higher backgrounds. 

We begin by discussing kinematics of the process. Consider a sequence of processes 

gg/qq^ X(q) ->• Vi(q 1 )V 2 {q 2 ), V x -> f(qu)f{qu), V 2 -> f(q 2 i)f(q 22 ), (1) 

that correspond to the production of a resonance X, followed by its decay to two vector bosons, followed by their 
decays to four fermions. The four-momenta of all particles are shown in parentheses. Momentum conservation implies 
Qi = <?ii + 1i 2 , q = <Zi + q 2 - We denote the invariant mass of the i-th gauge-boson by m| — qf and stress that it can 
differ from its mass m\. We assume that the particle X is produced on the mass shell, so that q 2 = {q\ + q 2 ) 2 = m 2 x . 
In what follows, we will refer to the heavier (lighter) of the two gauge bosons as V\ (V 2 ), mi > m 2 . 

As was already described in Ref. [20| . three invariant masses m VlV2 , mi, and m 2 , and six angles fully characterize 
the kinematics of the process in Eq. ([T]) in the rest frame of the resonance X. Five of these angles are illustrated in 
Fig. [TJ while the sixth angle defines the global rotation of an event in the plane transverse to the collision axis and, 
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FIG. 1: Illustration of a X particle production and decay in pp collision gg or qq — } X — > Vi (51)^2(52), VI — ► f (Q11) f (Q12) , 
V2 — » /(g2i)/(?22)- The three-momenta of the fermions (/) and antifermions (/), g n , q 12 , q 21 , and g 22 , are shown in their 
parent V% rest-frames, and the three- momenta of the Vi bosons, q 4 , are shown in the X rest-frame. For sign convention of the 
angles between planes see text. 



for this reason, it is not shown. We define these angles explicitly through the momenta of the leptons that are directly 
measurable experimentally. 

• The angles 8* G [0, tt] and $* G [— 7T, it] are defined through the unit vector of V\ direction, q x = 
(sin 9* cos sin#* sin cos 9*), in the rest frame of X. In this reference frame, the collision axis is aligned 
with the z-axis, h z = (0, 0, 1), taken as the direction of a colliding quark or one of the colliding gluons. Note, 
however, that the angle $* offset is arbitrary and it is not used in the final analysis. Also, when sequential decay 
of the vector bosons is not available, which is the case for X — > 77, only the angle 9* is accessible experimentally. 

• The angles $ £ [— ir,ir] and $1 € [— tt, tt] are the two azimuthal angles between the three planes constructed 
from the X decay products and the two V;-boson decay products in the X rest frame. The angle ^ € [— tt, tt] can 
be used in place of $1, it is defined as4' = $i + ( f>/2 and can be interpreted as the angle between the parton- 
scattering plane and the average between the two decay planes shown in Fig. [TJ These angles are explicitly 
defined as 

q 1 ■ (ni x n 2 ) -1 , ~ ~ \ 

$ = 1 x cos (-711 • n 2 ) , 

\q x ■ (m x 7i2)| 

q x ■ {hi x n sc ) 

*i = 1 ^-r: x cos ni • n sc , (2) 

\q 1 ■ (m x 7i sc )| 

where the normal vectors to the three planes are defined as 

<?11 X <?12 - 921 X <?22 , - ™z x q x 

"i = 1 1 , n 2 = 1 1 , and n sc = — . (3) 

\Qu x q 12 \ \<l 2 i x <l 22 \ |"*XQil 

In the above equations, q^^) is the three-momentum of a fermion (antifermion) in the decay of the Vi, and 
1i — In + 912 i s the Vi three-momentum, where all three-momenta are defined in the X rest frame. 



Finally, the angles 8\ and 9 2 e [0,7r] are defined as 



Ifc 9nl/ V I9ill9 2 i 



where all three- momenta are taken in the rest frame of Vi for the angle 9i . 
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The invariant masses of the two-fermion final states, the six angles defined above, and four-momentum of the initial 
partonic state exhaust the twelve degrees of freedom available to the four particles in the final state 3 . The initial state 
four-momentum defines the X invariant mass m VlV2 and the motion of the X system in the longitudinal (rapidity 
Y) and transverse (p T ) directions. Both Y and p T distributions depend on the production mechanism and therefore 
could help to further differentiate production models either for signal or background. However, these observables 
have little discrimination power between different signal hypotheses once production and decay channels are fixed 
and they introduce additional systematic uncertainties due to QCD effects. It is important to point out that the 
transverse momentum of the X particle introduces smearing in the determination of the production angles 9* and 
"J. The Collins-Soper frame [33| is designed to minimize the impact of the X transverse momentum on the angular 
measurements. However, the effect is expected to be small compared to statistical uncertainties for the luminosity 
expected in the 8 TeV run of the LHC and, for this reason, we do not study it in this paper. 

The full differential mass and angular distribution can be expressed using Eq. (|A1[) . where we can factorize the 
phase-space and propagator terms 



dT j (mi, rri2, cos 8*, 'J, cos 9\, cos #2, $) dTj(m±, m 2 , cos#*, cos#i, cos 62, $) 

oc 

drti\dra2 d cos 9* d^d cos 9\ d cos 82 d$ d cos 9* d^d cos 9\ d cos 92 d& 



xP(mi,m 2 ), (5) 



which are defined in Ref. [13[ as 

,J - - (mi-m 2 ) 



P(m 1 ,m 2 )= 1- (mi „ + ? m2) x 1- V '" L ^r" x — 2 „2iL m 2r2 x frr ,2 Z^rl , ■ ( 6 ) 
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v . 

(ml — m 2 , ) 2 + m 2 r 2 (m 2 — m 2 ) 2 + to 2 T 2 , 
After integration over the five angles, the differential mass distribution takes the form 



dTj 
dm\dm 2 



\A a p{m\,m%)\ x P(m 1 ,m 2 ). (7) 

a,/3=-,0,+ 



Below we discuss how to calculate Aq,^ (mi, 7712) for each spin and coupling hypothesis after a brief comment on 
the notation that we use throughout the paper. The polarization vectors of spin-one bosons are denoted by ef, we 
assume them to be transverse, q^i — 0. Fermion wave functions are conventional Dirac spinors. The spin-two X wave 
function is given by a symmetric traceless tensor t^„, transverse to its momentum ty, v q v = 0; its explicit form can be 
found in Ref. [2(j]- We will often use the notation /W^" = efqf — e^gf to denote the field strength tensor of a gauge 
boson with momentum qi and polarization vector e^. Assuming that momenta of the two bosons, Vi^, are along the 
z-axis qi } 2 = (£1,2, 0, 0, ±|q|), the polarization vectors read 

< a (0) = — (±|<?|,0,0,£i, 2 ), ei l (±) = e^(T) = ^=(0,Tl,-«,0). (8) 
1,2 V2 

The conjugate field strength tensor is defined as ffi} = 1/2 e^ va /3f^' a ^ — ^^va^q^- We use q — q± — 92 to denote 
difference of momenta of the two gauge bosons. 



A. Spin zero 



Suppose that the new boson is a spin-zero particle. The general scattering amplitude that describes the interaction 
of this boson with gauge bosons reads 

A(X H- V1V2) = UVmlelel+gWt^t m ^ , (9) 

where A denotes the scale where new physics could appear. We insert an explicit factor m\ in the amplitude to allow 
for a smooth massless limit consistent with generic requirements of gauge invariance which is relevant in case V = 7 
or g. 

It is instructive to discuss the connection between the amplitude in Eq. © and the concept of the effective La- 
grangian which is often used to discuss properties of the new boson. While the two approaches are related, the 



3 Throughout the paper, we take fermions in the final state to be massless. 



5 



amplitude A(X —> V1V2) provides a more general description of the properties of the new boson than any effective 
Lagrangian because the couplings are momentum-dependent form-factors that, for example, can have both real 
and imaginary parts. We do not expect this issue to be important for the new boson with a mass of 125 GeV, 
discovered at the LHC, but it may be essential for heavier resonances that may be discovered later, so we prefer to 
stick to this description. On the other hand, it is also true that effective Lagrangians lead to streamlined prediction 
for scattering amplitudes, since they provide an opportunity to order contributions of operators of different mass 
dimensions by their relevance, thereby reducing the number of terms that contribute to scattering amplitudes. Of 
course, given the scattering amplitude and assuming that form-factors are momentum-independent constants, the 
corresponding Lagrangian can always be constructed. For example, in case of Eq. (0), the following correspondence 
is valid 

(0) 2 (0) (0) 

V V M V 

g V) r W,^ f g)?!*r g f) Z , a Z^ [d p d a X] , g^f^f*^ ^ £ ~ g^XZ^Z^ (10) 

where v is the vacuum expectation value of the X held. Therefore, terms with g^ in A(X — > V1V2) are associated 

with dimension-three operators in the Lagrangian, terms with g^ 1 and with dimension-five, and terms with $3 
with dimension seven. As mentioned above, power-counting arguments suggest that lower-dimensional operators give 
larger contributions to the amplitude. 
We can re-write Eq. ([9]) as 

A(X -> ViV 2 ) = v^e^e™ (aig^m^ + a 2 q^q v + a 3 e MI , Q/3 qfq$ J . (11) 
The coefficients 01,2.3 are related to 5^,3,4 by 



ai = g 
where s is defined as 



,2 



s = qiq 2 =^ ^ 2 ~- (13) 

For a spin-zero resonance with couplings shown in Eq. (fTTj) . the three contributing helicity amplitudes are 

. m% ( n m 1 m 2 

A 00 = aiVl + ^ + a 2 5— x 



v 









mim 2 r- 






Ua 3 


2 v x 


V 






m x 




(a v 




mim 2 ,- 






-ia>3 


2 


V 






m x 



A ++ = 



A— = — ai-w 3 J— Vx) , (14) 



where x is defined as 



222 
m x — m{ — m 2 

2mym2 



(15) 



For a SM Higgs boson decaying to two massive vector bosons, ZZ or WW, the values of the couplings are g± = 1, 

and g^ = g^ — g± — 0. A small value of g 2 ^ ~ ©(chew) ~ 10~ 2 is generated in the SM by electroweak radiative 

corrections. The CP-violating constant gf is tiny in the SM since it appears only at the three-loop level. For the SM 

Higgs boson decays 77, Z7, or gg, only loop-induced couplings are possible so that g^' ^ while the other couplings 

are zero. However, allowing for beyond the SM scenarios, values of the gf^ need to be determined experimentally. 

For example, for a pseudoscalar Higgs boson one would expect g^ 1 ^ while the other g\^ = 0. It is also interesting 

to consider the model g^ ^ as an alternative to the SM scalar hypothesis, or a mixture of any of the above 
contributions. 
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B. Spin one 



For a spin-one resonance the amplitude depends on two independent terms 



A(X -> VtV 2 ) = b 1 [(e* iq )(e* 2 e x ) + (e^)(e^e x )] + b 2 e aiMU pe x e*^e* 2 ' v q 



(16) 



where e x is the polarization vector of particle X. The decay into two massless identical vector bosons is not allowed. 
The helicity amplitudes in the spin-one case corresponding to Eq. (1T6|) are the following 



.4, 



00 



(m| - mj) »_ 



-ib 



{ra\ — m 2 ) 



^4+o = bimi^fx + ih 



m 2 



A 0+ = -b\m 2 \fx-ib 2 



rrii 



A-o = b\rn\\fx — ib 2 



in- 

m 2 



(2 2 2\ / ''''i \ r\ 2 

m x — mj + m 2 ) — ^ — 1 ] + zrrijX 



( //r: I /;/ j - ///.] ) ( — | — 1 ) + 2m 2 X 



mr 



(ml _ TO 2 + m 2^ /'"h 2m\x 

* \ m 2 / 



mi 



/ 2 1 2 2\ 
(to x + m\ — m 2 ) 



1 ) + 2m\x 



(17) 



The model b\ = 7^ corresponds to a vector particle and b 2 = g^ 7^ to pseudovector particle, assuming 
parity-conserving interactions. Even though the spin-one hypothesis is rejected by the observation of X — > 77 decay, 
it is still interesting to consider the spin-one models in the decay to massive vector bosons. Indeed, there could be 
two nearby resonances at 125 GeV, one decaying to massive gauge bosons and the other to 77, and there have been 
models suggested [34| which predict the presence of two resonances. 



C. Spin two 



For a decay of a spin-two resonance to two vector bosons, including ZZ , WW, and 77, the scattering amplitude 
has the following general form 



A(X -> ViV 2 ) = A" 



rfin a 



(2)<T^\ f *(l) a f *(2) , 2 ( 9 (2) f I oJ2) ^1o: , , * v * a * a *v\ , (2) <P*T , * 

0,3 ' ^TiT 1 ^ V e l e 2 _ e l e 2 J + #7 V' 6 ! 



A 2 



A 2 



+ gs ^2 Z 11 '-' -'a/3 tm v I g9 ^2 e /^pc e l e 2 9 ' ^4 ^fiupcrl <1 k e l l</ e 2j+ e 2 \q t l)) 



(18) 



where t M „ is the AT wave function given by a symmetric traceless tensor [20] . This amplitude can be re- written as 



ci (giga)V + c 2 9^ut a 0q a q + c 3 g2M ^ liy t aP q a qP + 2c 4i q\ v q 2 t m + 2c 42 q2^q?t v 



A{X^V 1 V 2 )=K- 1 e*^e* 2 v 

+c 5 t a p- ! -^-€ t _ lupcr q^q 2 + c 6 t al3 q l3 e t _ iuap q p h ^ — {^a W aq p q J q v + t a v P oq p q a q^) 



(19) 



In case of massless bosons, like 77 or gg, the terms with m v in Eq. (1181) vanish. The coefficients ci_7 can be expressed 
through g{ 2 \ 10 



A 2 



= 2. 9 p) + 2 5 (^fl + ^i)fl + ^) + 2 5 (^, 
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C3 



C41 



C42 



A 2 ' 



Co 



= 9.(2) 




(20) 



We note that when constructing parametrizations of parity-odd amplitudes in Eq. ([15]). we should carefully exploit 
Schouten identities to remove mutually-dependent Lorentz structures. Such dependences lead to an interesting result 
- it turns out that a potentially contributing term t pv f£f va vanishes for traceless symmetric tensors, . 4 This 
cancellation implies that contributions due to t^ a q a e flupa e\'' 'e^ 'q a and V ia ei ia e^ pa q l 'e^ ^q a are related. Therefore, if 
we do not assume that the amplitude depends on only, we could have had two Lorentz structures contributing 
to the amplitude. However, because the Schouten identity connects these structures, we choose to keep only one of 
them in Eq. ([15]). 

We are now in a position to write down the helicity amplitudes for the spin-two case, using the parametrization 
shown in Eq. (|19p . For simplicity, we omit A in the following equations; the dependence on A can be restored on 
dimensional grounds. The amplitudes read 
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4 We are grateful to S. Palmer and M. Baumgart for useful discussions of this point. 
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A+- = A — l = m\% 



-C\X 



ij-miy 



ci 



(21) 



The minimal coupling scenario corresponds to the case g\ = g 5 7^ 0. However, when higher-dimension operators 
are considered, a broader range of options becomes available, analogous to g^ and in the spin-zero case. This 
variety of couplings corresponds to the complete set of vector boson V± and V% polarization states for the given mi and 

(2) 

m 2 - Among non-minimal couplings, the g\ term provides an interesting Lorentz structure with the field strength 
tensors of the two gauge bosons appearing similarly to the g% term in the spin-zero case. 

We note that, in principle, all couplings that we employ in the paper should be considered functions of kinematic 
invariants, e.g. g\ J ^ (m|, qf, Since the generic functional form of these couplings is unknown, accounting for 

dependences of g\ J ^ on q\ 2 introduces additional complications that are beyond the scope of this paper. Instead, 
we prefer to start the spin-parity determination program by treating couplings as constants to understand the "big 
picture" from data. Once this is accomplished, many further refinements and, in particular, kinematic dependences 
of the coupling constants, can be examined. We also note that we use the same parametrization for the amplitudes 
that describe decays X —> ZZ and X — > W + W~. In principle, since W's are not identical particles, the number of 
independent form factors required to describe the most general X — > W + W~ amplitude should be larger. We neglect 
this effect for the reasons explained above. Similarly, we point out that for spin-one and spin-two particles, the most 
general parametrization of the amplitude involves terms that depend on the difference of invariant masses of two vector 
bosons so that such terms vanish on the mass shell due to Bose symmetry. We do not include such terms in the present 
calculation and only employ those Lorentz structures in all amplitudes that have non-vanishing on-shell limit qf = 
q\ = m\. While this approximation is not parametric, we believe that the current parametrization already provides 
sufficient variety of Lorentz structures of couplings for hypothesis testing. More sophisticated parameterizations are 
only warranted if credible evidence shows that non-minimal couplings that we already introduced are insufficient to 
describe properties of the new particle. 



III. MONTE CARLO SIMULATION 

We have extended the simulation program [2(J |35[ to allow for various di-boson final states and to include the option 
of resonances decaying to off-shell gauge bosons. This program simulates the production and decay to two vector 
bosons of the spin-zero, spin-one, and spin-two resonances in hadron-hadron collisions, including all spin correlations. 
The processes gg/qq — > X — > ZZ and WW — > 4/, as well as gg/qq —> X —> 77, are implemented. It includes the 
general couplings of the X particle to gluons and quarks in production and to vector bosons in decay. The program 
can be interfaced to parton shower simulation (e.g. PYTHIA (36|) as well as full detector simulation through the Les 
Houches Event file format. 

As we discussed in Sec. HH in principle there is a large number of coupling constants to be determined. To illustrate 
the main idea of spin-parity determination, we pick several scenarios listed in Table [IJ Among them we include the SM 
Higgs boson spin-zero hypothesis (0+) and the graviton-like minimal coupling hypothesis for spin-two (2+). Other, 
more exotic, hypotheses are also considered. We note that for the spin-two scenarios, we assume that gluon fusion 



TABLE I: List of scenarios chosen for the analysis of the production and decay of an exotic X particle with quantum numbers 
J p . The subscripts m (minimal couplings) and h (couplings with higher-dimension operators) distinguish different scenarios, 
as discussed in the last column. The spin-zero and spin-one X production parameters do not affect the angular and mass 
distributions, and therefore are not specified. 



scenario 


X production 


X ->• VV decay 


comments 


0+ 


99 -» 


X 


g[ 0) / in Eq. © 


SM Higgs boson scalar 


K 


99 -> 


X 


g ( ° ] ^ in Eq. © 


scalar with higher-dimension operators 


<r 


99 -> 


X 


gf ] in Eq. © 


pseudo-scalar 


1+ 


qq -> 


X 


b 2 + in Eq. 


exotic pseudo-vector 


i" 




X 


61 ^ in Eq. 


exotic vector 


9+ 


g? ] + in 


Eq. JT8J 


gi 2) = gf 1 £ in Eq. (H8J 


graviton-like tensor with minimal couplings 


2+ 


g? ] + in 


Eq. (HHJ 


g{ 2) in Eq. (HHJ) 


tensor with higher-dimension operators 


2,7 


gW + in 


Eq. JTHJ 


gf ) ^ in Eq. JTHJl 


"pseudo-tensor" 
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FIG. 2: Distributions of some of the representative observables: mi in the X —> ZZ analysis (left), $ in the X — > WW 
analysis (middle), and cosO* in the X — > 77 analysis. Four signal hypotheses are shown: SM Higgs boson (red circles), 0~ 
(magenta squares), 2^ (blue triangles), 2^ (green diamonds), as defined in Table [I] Points show simulated events and lines show 
projections of analytical distributions. Here and throughout the paper, where only shapes of the distributions are illustrated 
and unless otherwise noted, units on the y axis are arbitrary. 

dominates the production mechanism, which is the case for the minimal coupling Kaluza-Klein graviton (2+ ) jl7j |. and 
this assumption may have an impact on the final results for the achievable significance of spin hypotheses separation. 
On the other hand, for the spin-zero scenarios, the production mechanism does not affect the angular and mass 
distributions. The chosen scenarios listed in Table U are similar to those considered in our earlier paper [201 ] . 

Distributions of some of the representative observables are shown in Fig. [5] for m x = 125 GeV. A complete set of 
distributions in the ZZ and WW final states is shown in Appendix [B] in Figs. [IT] [12] [13] Throughout the paper we 
consider yfs = 8 TeV proton-proton collisions and use the CTEQ6L1 parton distribution functions [371 ] . 

In the following we describe a simplified treatment of the detector effects which is not meant to reproduce exactly 
any of the LHC experiments, but still allows us to reliably understand feasibility of spin-parity studies at the LHC. 
We introduce smearing of the track momentum transverse to the collision axis, pt, and photon cluster energy. 
However, the exact resolution parameterization is not crucial as long as the overall signal-to-background separation 
is reproduced well. We mimic detector acceptance effects by cutting on geometric and kinematic parameters, such as 
Pt and pseudorapidity, r; = — lntan(#/2). Both leptons and photons are required to be in the effective acceptance 
range |?7| < 2.5. 

The main backgrounds in the X — > ZZ, WW, and 77 analyses are the continuum di-boson production, includ- 
ing Z7* for ZZ These are modeled with POWHEG ^ (ZZ) and MadGraph [H (WW, 77). Additional 
contributions of backgrounds with fake vector boson reconstruction requires special treatment. However, their con- 
tributions are smaller and observable distributions are similar to the VV background, so their contributions can be 
effectively accounted for by rescaling the di-boson background rate to match total background rates observed by the 
LHC experiments. 

IV. ANALYSIS METHODS 

In this Section, we illustrate the application of the matrix element analysis formalism to distinguishing different 
spin-parity hypotheses for the observed boson near 125 GeV. We illustrate this with the seven scenarios defined in 
Table U and comment on future direction of the measurements. 

In Ref. [20( we pointed out that the ultimate goal of the analysis should be the experimental determination of all 
helicity amplitudes that involve X and two gauge bosons. The techniques discussed here and in Ref. [20j are ideally 
suited for such measurements since parameters in the angular and mass distributions become fit parameters in analysis 
of data. However, such multi-parameter fits require large samples of the signal events which are not yet available. 
Therefore, in our opinion, the first step in understanding the spin-parity of the resonance should be distinguishing 
between different hypotheses. For such a goal, a simplified, but still optimal, analysis approach can be developed that 
employs just two observables. A simple extension of this analysis, which naturally arises if we assume, for example, 
that the observed resonance is a mixed spin-parity state, is to fit for ratios of couplings. Ultimately, this approach 
will lead to a complete multi-dimensional fit of all coupling parameters using a complete set of kinematic observables. 
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TABLE II: Expected separation significance (Gaussian a) between the SM Higgs boson scenario (0„) and various hypotheses 
defined in Table [T] Expectations are given for the scenario of a 5.0 a signal-to-background separation observed in the search for 
the SM Higgs boson in each channel, and therefore interpretation in terms of integrated luminosity and pp collision energy at 
the LHC may differ significantly between the three channels X — » ZZ, WW , and 77. 



scenario 


X ^ ZZ 


X -> WW 


X — > 77 


0m vs background 


5.0 


5.0 


5.0 


0+ vs 0+ 


1.7 


1.1 


0.0 


0+ vs 0- 


2.9 


1.2 


0.0 


0+ vs 1+ 


1.9 


2.0 




0+ vs 1- 


2.6 


3.2 




0+ vs 2+ 


1.5 


2.8 


2.4 


0+ vs 2+ 


~5 


1.1 


3.1 


0+ vs 1~ 


~5 


2.5 


3.1 



Going back to the two-dimensional fit, we note that one of the two observables is related to the resonance mass as 
it typically has the largest discriminating power against the background. This observable depends on the final state; 
for example, it is the four-lepton invariant mass ma in the X — > ZZ — > A£ analysis, the transverse mass tut [J 0] in 
the X — > WW — > 2£2v analysis, or the two-photon invariant mass m 77 in the X — > 77 analysis. 

The second observable combines other kinematic information that is available, and it is designed to distinguish 
between different signal spin-parity hypotheses in the optimal way. In the X — > ZZ — > 4£ analysis we build the 
kinematic discriminant, defined in the MELA approach adopted by the CMS experiment @, [25j], which combines the 
five angular and two mass observables in the optimal way. In the X — > WW — > 2£2v analysis, the complete matrix 
element information cannot be exploited because of the neutrinos in the final state. Therefore, we adopt a simplified 
approach by picking one observable that is most sensitive to the spin-parity of X. We found this observable to be the 
di-lepton invariant mass m« while the opening angle between the two leptons in the transverse plain provides less 
sensitivity. Finally, in the X — > 77 analysis, the only available observable is cos 8* since there is no further sequential 
decay chain involved. 

We note that it is not our goal in this paper to demonstrate how the analysis should be optimized for the signal- 
background separation. Doing so requires simulation of detector performance and of all background processes 0, Q ■ 
Instead, we assume an excess of signal events over background in each of the three channels, X — > ZZ, WW, and 77, 
and calculate the achievable level of separation power between different signal spin-parity hypotheses. While precise 
prediction of spin-parity separation significances also requires detailed simulation, as long as the phase-space of the 
discriminating observables is well-modeled, such predictions are less sensitive to details of the analysis once a given 
signal significance is observed. 

We present results for the expected separation significance between the SM Higgs boson scenario (0+ ) and various 
J p and coupling hypotheses defined in Table U for a fixed hypothesis of a signal excess, which we take to be 5a for 
the SM Higgs- like resonance. The performance quoted in Tabic [TTI follows from the studies presented in the following 
subsections and can be interpreted in terms of integrated luminosity and pp collision energy at the LHC for each of the 
three channels X — !• ZZ, WW, and 77. We observe that a simple rule of scaling with luminosity L, significance ~ \fh, 
is a very good approximation in these studies as long as the uncertainties are dominated by statistical errors. 

We use an extended maximum-likelihood fit [2(| to extract simultaneously the signal and background yields. The 
likelihood is defined as 

C k = cxp (-n s i g - nbkg) J} (n sig x V^ g (xi; a; j3) + n hkg x Vbk g (xi, (3)) , (22) 

i 

where n S i g is the number of signal events, nbkg is the number of background events and Vixi; a; (3) is the probability 
density function for background or signal for different spin hypotheses, k. Each event candidate, i, is characterized 
by a set of two observables Xi = (m,D). The signal coupling parameters are collectively denoted by a, and the 
remaining parameters by (3. The correlated (m, D) distribution is parameterized with a binned histogram (template) 
using simulation. The likelihood Ck in Eq. (|2"2"]l is evaluated independently for each spin hypothesis k. Two sets of 
pseudo-experiments are generated, each with same average number of signal events of a particular type embedded 
into the expected background. 

Examples of distributions of — 21n(£i/£2) are shown for a large number of generated experiments in Fig. [3l where 
one of the signal types is chosen to be the SM Higgs boson. The probability for an alternative signal to produce a 
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FIG. 3: Distributions of — 21n(£i/£2) with the likelihood Ck evaluated for two models and shown for a large number of 
generated experiments in the analysis of X — > ZZ (left), WW (middle), and 77 (right). The models shown are the SM Higgs 
boson + (red solid points) and the pseudoscalar 0~ for ZZ or the graviton-like 2„ for WW and 77 (blue open points). The 
scenarios correspond to those shown in Table HTl 



value of — 2ln(Ci/ C2) below the median value for the SM Higgs boson hypothesis is taken as the one-sided Gaussian 
probability and interpreted as the number of Gaussian standard deviations, S. The value of S corresponds to an 
effective separation between the two distributions in the symmetric case, or equivalently to the expected separation 
between the two hypotheses. However, a certain amount of asymmetry between the distributions is possible, as we note 
in some cases below, and the expected significance of separating type 2 signal from type 1 may differ from separation 
of type 1 from type 2. An approximate average of the two values could be obtained from the point beyond which the 
right-side tail of the left histogram and the left-side tail of the right histogram have equal areas (corresponding to 
5/2). We choose to quote the first of the three values as more relevant for separation of alternative hypotheses from 
the SM Higgs boson. A similar technique can be employed for the significance calculation of the signal excess over 
background. Below we discuss details of the analysis methods that are particular to each channel. 



A. X — >■ ZZ 



In the X — > ZZ —> 4£ channel, the dominant background is the continuous production of Z-f*/ZZ 5 . The ZZ 
production cross section is comparable to that of the SM Higgs boson in the four-lepton invariant mass window 
comparable to detector resolution. Loose selection requirements are applied to simulated signal and qq — > ZZ 
background events to model detector effects of CMS and ATLAS. For lepton track transverse momentum, we apply 
Gaussian random smearing with an rms Apr — 0.014 x pr (GeV) for 90% of the core of the distribution and a wider 
smearing for the 10% tail. Leptons are required to have pseudo-rapidity in the range |?7| < 2.4 and pr greater than 
7 GeV. In addition, leptons with the highest and next-to-highest transverse momentum are required to also have 
Pt > 20 and 10 GeV, respectively. To reject the non-ZZ background, the invariant masses of the di-lepton pairs are 
required to satisfy 50 < toi < 120 GeV and 12 < m-i < 120 GeV, where mi > m^. The overall ZZ rate is then 
scaled to be consistent with the total background observed in LHC experiments, including Drell-Yan and top events 
with jets faking leptons. We do not attempt to model this instrumental background shape and implicitly assume that 
shapes are well-modeled by qq — > ZZ events. Only events with 110 < < 160 GeV are considered in the final 
analysis. The number of signal events after all selections is 0.8 events/fb _1 , while the number of background events 
is 1.9 events/fb" 1 . Using only the m<u shape of signal and background, we find an expected significance of 3.3 a with 
10 fb _1 of data, comparable to that observed at the LHC 6 . 

As was pointed out earlier in Ref. (20l . [23j . using full kinematic information in the X — > ZZ channel improves 
signal-to-background separation by about 20% compared to a one-dimensional analysis of the invariant mass m«. 
This has been exploited by the CMS experiment in the discovery of the new boson One can perform either 



5 We will collectively refer to these two processes as ZZ in what follows. 

6 We disregard the difference between the 7 and 8 TeV collision energies of LHC for simplicity. 
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a multidimensional fit or create a kinematic discriminant (MELA) 0, which is constructed from the ratio of 
probabilities for signal and background hypotheses 



Dbk B = 



V s - ig {mu] mi, m 2 ,0) 



(23) 



Here V s i g and "Pbkg are the probabilities, as a function of masses nij and angular observables ft for a given value of 
invariant mass as defined in Eq. ([5]), for the SM Higgs boson signal and ZZ background, respectively. Although 
analytic computation of the matrix element for continuum ZZ production [23| exists, it neglects the Z^j* process and 
therefore it can not yet be applied to the region below mu ~ 180 GeV. Instead, we use a large sample of POWHEG 
simulated events to fill a multi-dimensional histogram (template), where the most important correlations between up 
to three observables are taken into account. 

The above approach to background rejection is illustrated in Fig. 21 where we plot mu and -Dbkg, which are 
mostly uncorrelated in the small ma region considered. As shown in Fig. 21 for a wide range of different signal 
spin-parity hypotheses, the £>bkg distributions do not differ considerably. However, all signal distributions of -Dbkg 
differ considerably from background. We confirm that significance of the signal observation in the two-dimensional 
analysis of (mui -Dbkg) increases by more than 20% compared to a one-dimensional analysis of the invariant mass 
run- To simplify the fitting model, in the rest of the paper we will not use the additional background suppression 
power of the Dbkg observable but we note that the effective significance can be increased by either including Dbkg in 
the multivariate fit, or, equivalently, using information contained in Dbkg in the fit. 

Separation between different spin-parity scenarios of the observed resonance can also be obtained using kinematic 
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FIG. 4: Distributions of rriAe (left), -Dbkg (middle), and D - (right) in the X — > ZZ analysis for the non-resonant ZZ 
background (black solid circles), and four signal hypotheses: SM Higgs boson (red open circles), 0~ (magenta squares), 2^ 
(blue triangles), and 2^ (green diamonds). Not all signal hypotheses are shown on all plots. The mass range 120 < ni4£ < 130 
GeV is shown in the D distributions. 
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FIG. 5: Template distributions of D - vs ma in the X — > ZZ analysis for the SM Higgs boson (left), pseudoscalar resonance 
(middle), and non-resonant ZZ background (right). 
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information. We can re-write Eq. (|23j) as 



D, 



j ?2(w4<;mi,m2, n) 
7>i (m 4 £;mi,m 2 ,r2) 



(24) 



where Pi and 7^ are the probabilities as defined in Eq. ([5]), for two different hypotheses of spin-parity and tensor 
structure of interactions of the signal resonance. Equation (|24p is indeed the optimal way to combine all relevant 
kinematic information into a single observable Djp for separating the SM Higgs boson scenario from other 
hypotheses without loss of information. 

As an example, in the right plot of Fig. U we show a kinematic discriminant D - optimized for separating 0+ and 
0~ signal hypotheses. D - is calculated in Eq. (|24|) with V2 taken as the probability density for - . Since mu is still 
the most powerful observable to discriminate any type of signal from the background, we perform a two-dimensional 
fit of (mu, D -). The probability densities for signal and background are parameterized as two-dimensional template 
histograms using simulation, as shown in Fig. [S] 

In the scenario described above, with an expected SM Higgs boson signal significance of 3.3 a with 10 fb -1 of data, 
we estimate an average separation of 1.9 a between the 0+ and CP signal hypotheses. Equivalently, assuming that the 
integrated luminosity is high enough to ensure 5cr signal-to-background separation, the average expected separation 
of 0+ and 0~ is 2.9 a, see Fig. [3] This and other results for several other signal hyp otheses are shown in Table ITT1 We 
find the 0+ and 0~ separation results consistent with those predicted by CMS [40(, taking into account the assumed 
signal significance (expected vs observed). 

The separation power depends on information contained in kinematic distributions; we show illustrative examples 
in Figs. [U QTJ and [T2J For example, separation of the SM Higgs boson hypothesis from 0~ is better than from 2+ 
since a number of mass and angular distributions are more distinct. Also, we note that both 2^ and 2^" are even more 
different from 0^ than any other hypothesis considered. One of the kinematic distributions that shows important 
differences between the SM Higgs boson and the pseudoscalar, as well as between the SM Higgs boson and 2^ or 2^, 
is the low-mass tail of mi distribution, see the left plot in Fig. [5] As shown in Ref. [4(|, there is a rather large fraction 
of the X — > ZZ events in CMS with both Z's off-shell. If this feature persists in data, it may reveal contributions of 
more exotic couplings shown in Eqs. © and (|18p . We note that a particular feature which may enhance the mi tail 
considerably, as shown for the 2^~ hypothesis in Fig. [21 is the presence of a large power of the parameter x in Eq. (|2ip . 



(2) 

In turn, the appearance of this parameter is related to the terms C2 and C3 in ^4oo which are sensitive to gl coupling. 

We comment about an interesting feature in Table HU On general grounds one can expect that significance of 
hypotheses separation between two types of signal to be smaller than the observation significance of the signal. 
However, in a situation when the kinematic discriminant itself provides substantial background rejection power, 
significance of the observation of the alternative signal may become higher than that for the SM Higgs-like resonance. 
This phenomenon occurs in the study of 2^ and 2^" hypotheses, where the m\ mass distribution becomes a particularly 
powerful discriminating observable. As a result, for the corresponding signal types, D 2 + and D 2 - may become even 

^ h h 

stronger background rejection observables than mu- We do not see this in Fig. @] because the SM Higgs boson 
hypothesis is used for the computation of -Dbkg, but an alternative signal hypothesis could have been considered as 
well. The mi-distribution then also leads to a very strong signal hypothesis separation, approaching the values of 5 a 
for SM Higgs boson vs 2\ and 2^ in Table HO 

We note that analysis of the spin-parity hypotheses should not be limited to just discrete hypothesis testing. In 
Ref. [20j we showed how a continuous spectrum of parameters can be obtained from a multidimensional fit. As an 
intermediate step, one could consider determination of the fraction of a certain component in a mixed state. For 
the spin-zero particle, this can be modeled by non- vanishing g{ ^ and g{ ^ couplings in Eq. ([5]). Wc note that, in 
this scenario, there is an interference term in the amplitude which is correctly described by our simulation. Results 
presented in Table [TT] can be used to illustrate the typical precision on the fraction by dividing the full range between 
the two extreme hypotheses by the number of standard deviations between them. For example, this implies that by 
the end of the 8 TeV run, the LHC experiments may be able to constrain an approximately 50% admixture of the 
CP-violating amplitude at 95% confidence level. 



B. X WW 



Compared to the ZZ final state, the X — > WW — > 2£2v channel is expected to have a larger rate due to a larger 
branching fraction of WW — > 2l2v 1 provided that decay rates X — > ZZ and X — > WW are comparable. However the 
analysis suffers from large backgrounds and the fact that neutrino momenta cannot be reconstructed. In place of the 
four-lepton invariant mass, the transverse mass defined as tot = (2p^E^ lss (l — cos Acj) u _ E n^ B!1 )) 1 ^ 2 , where A(f> u _ E n,i BS 

is the angle between the direction of the di-lepton pair and the missing energy _E™ 1SS vector in the transverse plane, 
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FIG. 6: Distributions of mr (left), mu (middle), and A<f)u (right) in the X — >• WW analysis for the non-resonant WW 
background (black solid circles), SM Higgs boson (red open circles), and a spin- two resonance in the 2+ model (blue triangles). 
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FIG. 7: Template distribution of mu vs uit in the X — >• WW analysis for the SM Higgs boson (left), spin-two resonance 
hypothesis 2+ (middle), and non-resonant WW background (right). 



is exploited to disentangle signal from background [HQ- ^ n our simplified study _E™ 1SS is calculated from the 2u 
momentum. The one-jet e/j, and zero-jet same-flavor categories only contribute to the signal sensitivity at the 10% 
level because of larger backgrounds from top-quark decays and Drell-Yan production, respectively 0, • Therefore, 
we only select events with different lepton flavors (e/i) and little jet activity to enhance signal-to-background ratio, so 
all events with jets with transverse energy greater than 30 GeV are rejected. In this category, the main background 
comes from the non-resonant WW production 0, Q . To further reject the reducible backgrounds such as Drell-Yan 
and W + jets/7 processes, we require pr > 20(10) GeV for the leading (sub-leading) lepton, p^ > 30 GeV for the 
di-lepton system, E™ iss > 20 GeV, 60 < m T < 130 GeV, and 10 < m« < 90 GeV. 

We estimate the expected number of signal and WW background events after this selection by extrapolating the 
expected yields in the signal regions used in Ref. Q to the signal region denned above using simulation. The estimated 
number of SM Higgs boson events is 13 per fb _1 . The number of non-resonant WW background events is estimated 
to be 104 per fb _1 . We also assume that continuum WW production gives two-thirds of the total background 
and that kinematic distributions of the non-IUIU backgrounds are the same as the ones of the WW background. 
We cross-check this estimation using the signal region used in Ref. [l[ and find consistent results. To extract the 
expected significance for separating different signal hypotheses S, we construct a two-dimensional template based on 
two observables (m,D) = (%,m«); this is illustrated in Figs. [B] and [7] We have also considered other observables, 
such as the azimuthal angle A<^« between the two leptons and found smaller separation compared to the case when 
ma is used. On the other hand, since there is large correlation between A</>« and mu, using three observables in the 
fit is not expected to increase the significance of the separation much. 

Using this simplified background model, we estimate the expected significance for distinguishing the SM Higgs 
boson hypothesis from the background with 10 fb _1 using either the single observable mu or the two observables bit 
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and mu. The former approach gives 2.6 a separation from the background, similar to results of the LHC |l|, |2j, while 
the latter gives 3.5 a which is an improvement of 35%. We follow the procedure outlined for the X — > ZZ analysis 
above and present the results in Fig. [3] and Table ITll We find good separation between the SM Higgs boson and the 2+ 
hypotheses in particular, where this channel may have an advantage over the X — > ZZ channel. The reason for better 
performance of the WW channel for 2+ separation is the larger value of parameter A f defined in the Appendix [K\ 
which enters the angular distributions in Eq. (|A1|) . As a consequence, there are larger azimuthal angular variations 
which are illustrated in Figs.[31[T21 andlT3"l 

C. X ->• 77 

In the inclusive X — > 77 decay analysis, all information about the couplings is contained in the cos 9* distribution. 
The distribution is flat for a spin-zero resonance, while for a spin-two it is a normalized second degree polynomial in 
cos 2 9* which requires two independent parameters. Non-zero values of either parameter would be an unambiguous 
sign of a spin-two (or in principle higher spin) resonance. However, relating these coefficients to general couplings will 
have many ambiguities which are not generally present in the ZZ and WW channels. Indeed, the spin-two X — > 
angular distribution reads 




m yy [GeV] cose* 



FIG. 8: Distributions of m 77 (left) and cos#* (right) in the X — > 77 analysis for the non-resonant 77 background (black 
solid circles), SM Higgs boson (red open circles), the spin- two resonance in the 2^ model (blue triangles) and 2^ or 2^ models 
(green diamonds). The mass range 120 < m 77 < 130 GeV is shown in the cos 6* plot for background. 




FIG. 9: Template distributions of cos#* vs m 77 in the X — > 77 analysis for the background (left), SM Higgs boson (middle), 
and spin-two resonance hypothesis 2+ (right). 
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16 ^ ( f. 7 " 2 ^ 77) = (2 - 2f zl + f z2 ) - 6(2 - Af zl - f z2 ) cos 2 9* + 3(6 - 10/,! - 5f z2 ) cos 4 9* 

1 (X COS U 

+/+_ {(2 + 2/,i - 7/ z2 ) + 6(2 - 6/,i + /, 2 ) cos 2 0* - 5(6 - 10/ 2l - 5/ z2 ) cos 4 0*} 

(xl + Ax cos 2 9* + B x cos 4 9* , (25) 

where /++, / , and / |_ = /_| are fractions of transverse amplitudes in the decay, and f z \ and f z2 are polarization 

fractions in production, see Appendix A of Ref. [20] for more details. The special case of the minimal coupling 

in both production and decay corresponds to f z \ + f z2 — 1 and /_| = / |_ = 1/2. In this case, one obtains 

(1 + 6cos 2 6>* + cos 4 0*) for the gg production mechanism with f z2 = 1 and (1 — cos 4 9*) for the qq production 
mechanism with f z \ — 1. The ideal distributions in Eq. (|2"5)) are shown together with generated events in Fig. [5] 
These distributions are identical for the 2^ and 2^ hypotheses. 

For illustration purposes, we proceed with the discussion of a simplified analysis. The acceptance thresholds, chosen 
to be similar to those used in LHC analyses, are E\ > m 77 /3 and Ej, > m 77 /4 for the first and second photons, 
respectively. We apply ^-dependent Gaussian random smearing to photon cluster energy which varies between between 
1% in the central pseudorapidity region and 6% in the forward region. Photons must be inside the calorimeter 
acceptance |f?| < 2.5 and outside the crack region 1.44 < \r)\ < 1.57, similarly to the CMS experiment. Considering 
the mass window 110 < m 77 < 140 GeV, we estimate the expected number of signal and background events based on 
Refs. [l|,[l] to be 22 and 3515 per fb _1 , respectively. To extract the expected separation significance between different 
signal spin-parity hypotheses S, we construct a two-dimensional template based on (m, D) = (m 77 ,cos#*) shown in 
Figs. [5] and [SI The loss of events at large values of |cos0*| limits the precision of polarization measurements and is 
due to px and rj selection requirements. Similar effects appear in the analysis of the Drell-Yan process, as discussed 
for example in Ref. [4l|. We rely on the shapes of the distributions after the above kinematic selection, and the 
normalization is taken from data. Using two-dimensional (m 77 ,cos#*) templates, we obtain 2.7 a significance with 
10 fb _1 , which is similar to the LHC results expected for the SM Higgs boson [H, 0). 

We find good separation between the SM Higgs boson hypothesis and the spin-two models considered, as can be 
seen in Fig. [3] and in Table [TTl However, since just one angle is available in the analysis, the separation power may be 
weak or absent for other models where the cos 9* distribution is close to flat. 



V. SUMMARY AND CONCLUSION 



We have described a framework to determine the spin, parity, and general tensor structure of interactions of the 
new boson observed at the LHC. We consider a variety of Lorentz structures for spin-parity hypothesis testing that go 
beyond the minimal couplings expected for the SM Higgs boson or the graviton-like interactions of a spin-two boson. 
The full analytical calculation of angular and mass dependence of the decay amplitude X — > V*V* allows the most 
general analysis of a resonance with any integer spin J. A Monte Carlo simulation of the process pp — >• X — > V*V* , 
with V = Z, W, and 7, with off-shell electroweak gauge bosons, all spin correlations, and general couplings enables 
experimental investigation of the properties of the new resonance. Both the analytic formulas and the event generator 
are publicly available, see Ref. (35[. 

We have illustrated how the spin and parity of the new boson can be tested in the processes pp — > X — > ZZ, 
WW, and 77, using simplified simulation of the background and of the detector effects at the LHC experiments. We 
have presented the expected significance of spin-parity hypothesis separation for several scenarios in Table UH where 
it is assumed that the 5a signal-to-background separation is achieved in each channel. The linearity of the relation 
between the signal-to-background significance and the spin-parity signal hypothesis separation significance allows us 
to extrapolate expectations to different luminosity scenarios, as shown in Fig. [TU] for 0~ and 2+ models. We rely 
on the expected signal-to-background significances reported by the LHC experiments for the integrated luminosity 
of about 10 fb -1 , which we take as 3.8, 2.4, and 2.8 a in the X — > ZZ, WW, and 77 channels, respectively @. In 
Table Hill we show examples of hypothesis separation expectations, per each LHC experiment, by the end of the 8 TeV 
LHC run, assuming 35 fb" 1 of integrated luminosity. 

We would also like to comment on some other potential final states in the decay of the new boson, such as Z"f and 
fermion-antifermion final states. Since no significant excess of events in these final states has been observed, we leave 
detailed discussion of these final states to later work. However, the techniques discussed in this paper are applicable 
to them as well. For example, analysis of associated production qq —> Z* — >• ZX —> (£~£ + )(bb) or (£~£ + )(t~t + ), 
and similarly qq' — > W* — > WX, would follow the same formalism as discussed above. In the above processes, the 
angular distributions of decay products should allow discrimination between the spin and coupling hypotheses for 
both XVV and Xff. For a spin-zero X decay, the fermion angular distributions are flat; for a spin-one X, the 
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TABLE III: Expected separation significance S (Gaussian a) between the SM Higgs boson scenario (0^) and 0~ or 2+ 
hypotheses in the analyzed channels and combined, for the scenario corresponding approximately to 35 fb~ of integrated 
luminosity at one LHC experiment. 



scenario 


X -> zz 


X -> WW 


X -¥ 77 


combined 


0^ vs background 


7.1 


4.5 


5.2 


9.9 


0+ vs 0" 


4.1 


1.1 


0.0 


4.2 


0+ vs 2+ 


2.2 


2.5 


2.5 


4.2 




3 4 5 6 
0+ vs bkg [a] 



3 4 5 6 

0+ vs bkg [a] 



FIG. 10: Expected hypotheses separation significance vs signal observation significance for the SM Higgs boson vs 0~ (left) 
and 2+ (right) hypotheses. Points show two luminosity scenarios tested with generated experiments and expectations are 
extrapolated linearly to other significance scenarios. Dashed lines indicate what might be expected with 35 fb _1 of data at one 
LHC experiment. 



angular distributions are similar to those in Eq. (|A1[) with J = 1; and for a spin- two X, the angular distributions 
can be obtained in a similar manner. The angular distributions for general couplings in the decay X — > ff can be 
obtained from Rcf. 20]. 

In view of the importance of the discovery of the new boson for particle physics, it is important to confront all 
theoretical assumptions about its properties against experimental facts. The results presented in this paper point 
towards a realistic possibility that by the end of the 8 TeV run of the LHC, extreme hypotheses about spin and 
parity of the new boson can be experimentally excluded. However, it will be much harder to exclude contributions of 
anomalous couplings of the Higgs bosons to gauge bosons if they are smaller than ten percent of the SM couplings. 
For that, a significantly larger dataset will be required and multivariate fitting techniques [20] will help to achieve this 
goal. Nonetheless, if the nature of the new boson discovered at the LHC is exotic, there is a good chance to determine 
this already in the coming year. 
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Appendix A: Analytical angular distributions 



In this appendix we present the general angular distribution in the production and decay of a particle X, with any 
integer spin J, in parton collisions ab -t X -t Vi(qi)V 2 (q2), V\ ~> /(«ii)/(«i2). V 2 -+ /feO/Cfe), as derived in 
Ref. [20| and generalized here to remove the constraint between the A a p and Ap a amplitudes. Helicity amplitudes 
A a p depend on the vector boson resonance masses mi and m.2, as described in Eqs. (Tl4|) . (fT7|) . and ([21~j) . and related 
formulas incorporating the couplings. We work in the rest frame of the resonance X and all angles that we use below 
are defined in Section II. 

The amplitudes A a p are, in general, complex and the angular distribution is parameterized by the magnitude of 
the amplitude \A a p\ and the phase </> a p — arg(A a p / Aqq) . The angular distribution is 

Afj dT J (m 1 ,m 2 , cos 0* cos 6* 1; cos 2 , $) _ 



d cos 6* <N>d cos 0i d cos 2 d$ 
N ) x U |A 00 | 2 sin 2 0i sin 2 2 



+ \A ++ \ 2 (l + 2A fl cos 6»i +cos 2 6 X ) (l + 2A h cos0 2 + cos 2 2 ) 
+ |A__| 2 (l - 2% cos 6»i + cos 2 0i) (1 - 2A h cos0 2 + cos 2 2 ) 
+ ^^oollAf+KA/i + cos 0i ) sin0i(A /2 + cos0 2 ) sin0 2 cos($ + <j> ++ ) 
+ 4\A ao \\A | ( A fl - cos 0i ) sin X ( A h - cos 2 ) sin 2 cos($ - 4> ) 

+ 2|A ++ ||A | sin 2 0i sin 2 2 cos(2$ - cp + 4> ++ ) 

*) x [2|A +o | 2 (l + 2A /l cos0i+cos 2 0i)sin 2 2 
+ 2|A -| 2 sin 2 0i (1 - 2A h cos 2 + COS 2 ) 
+ 2|A_ | 2 (1 - 2% cos 0i + cos 2 0i ) sin 2 2 
+ 2\A 0+ \ 2 sin 2 0i (1 + 2A h cos0 2 + cos 2 2 ) 

+ 4\A +0 \\A -\(A fl +cos0i)sin0i(A /2 -cos0 2 )sin0 2 cos($ + </> +0 - 4>a-) 
+ 4| J 4 0+ ||A_ |(A /l -cos0i)sin0i(A /2 4-cos0 2 )sin0 2 cos($ + <f> a+ - 4>-o) 

+Fl_ 1 (6*) x 4|A +0 ||A 0+ |(A /l +cos0i)sin0i(A /2 + cos0 2 ) sin0 2 cos(2* - +o + </> 0+ ) 

+ 4|A +0 ||A_o| sin 2 0i sin 2 2 cos(2# - $ - +o + 0_o) 
+ 41^0-11^0+ 1 sin 2 0i sin 2 2 cos(2W + $ - O _ + O+ ) 

+ 4| A _ | |A_ 1 ( Afi - cos #i ) sin 0i /2 ~ cos 0a) sin 2 cos(2* 
f ) x |~|A + _| 2 (1 + 2A h cos 0i +cos 2 0i)(l - 2A/ 2 cos0 2 + cos 2 2 ) 



+^ 2 J 2 ( 



|A_+| 2 (1 - 2A h cos 0i + cos 2 0i)(l + 2A/ 2 cos0 2 + cos 2 2 ) 



+F 2 J _ 2 (9*) x 2|A + _P_+|sin 2 0isin 2 2 cos(4* - 



-) 



+F o J i(0*) x ^V^ooP+oK^ +cos0i)sin0isin 2 2 cos(*-$/2-0 +o ) 

+ 4n/2|A 00 ||v1o-| sin 2 6 x {A h ~ cos0 2 ) sin0 2 cos(* + $/2 - <fo-) 

+ 2V2|A ||Afo| sin 2 0i(A /2 - cos 2 ) sin 2 cos(-* + 3$/2 + </> +0 - 

+ 2V2\A— \\A - \{A h - cos 0i) sin 6 X (1 - 2A /2 cos 2 + cos 2 2 ) cos(-# + $/2 + <fo- 
+ 2\/2|A+ + ||A +0 |(l + 2A/j cos0i + cos 2 0i)(A /2 + cos0 2 )sin0 2 cos(* + $/2 - </> +0 
+ 2v / 2|A ++ ||A _|(A /l + cos 0i) sin 0i sin 2 2 cos(* + 3$/2 - + <?(>++) 

+^-l(0*) x [4V2|AooPo + |sin 2 0i(A /2 +cos0 2 )sin0 2 cos(* + $/2 + </) + o) 
+ 4V2| Aio 1 1 A- \ (A h - cos 0i ) sin d x sin 2 2 cos(* - $/2 + 0_ o ) 
+ 2v / 2|A__||A 0+ |(A /l - cos 0i ) sin 0i sin 2 2 cos(* + 3$/2 + O + - 0—) 
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-F, 



0,2V 



+ 2\/2|yl P_ |(l - 2A h cos 9\ H- cos ^1) (-A/2 — cos 62) sin #2 cos(^ -h $/2 -h 0_ — ) 

+ 2V2|A ++ Po+|(% +cos0i)sin0i(l + 2A/ 2 cos0 2 + cos 2 2 ) cos(* - $/2 + </> + - 
_o| sin 0i(A/ 2 + cos0 2 ) sin0 2 cos(-* + 3$/2 - + 

V) x 4|v4 o||A + _|(A /l +cos0i)sin0i(A /2 - cos 2 ) sin 2 cos(2* - <j> + J) 

+ 2|A 1 1 A+_ I sin 2 0i (1 - 2A h cos 2 + cos 2 2 ) cos(2* - $ + - <?!>+_ ) 

+ 2| I A+_ |(1 + 2% cos 0i + cos 2 X ) sin 2 2 cos(2* + $ + ++ - 0+_ ) 

+ F o 7 -2(^*) x 4 IA)o||^-+|(% -cos0i)sin0i(A /2 + cos 2 ) sin 2 cos(2* + </>_+) 

+ 2|A 1 1 A_ + |(l - 2% cos 0i + cos 2 0i) sin 2 2 cos(2* + $ - + 0_+) 

+ 2 1 A ++ 1 1 A_+ 1 sin 2 X (1 + 2A/ 2 cos 2 + cos 2 2 ) cos(2# - $ - cj> ++ + 0_+ ) 

*) x 2- v / 2|A K) ||A l __|(l + 2A h cos0i + cos 0i)(^4/ 2 — cos0 2 ) sin0 2 cos^ + $/2 -+ 

+ 2\/2\A Q ^\\A + -\(A fl + cos0i) sin 0i (1 - 2A h cos0 2 + cos 2 2 )cos(* - $/2 4 

- 2%/2|^l 0+ P_ + |(A /l - cos0i) sin 0i (1 + 2% cos0 2 + cos 2 2 ) cos(-* + $/2 

- 2- v / 2|A-o||-4-+|(l - 2% cos 0i + cos 2 0i)(A /2 + cos0 2 ) sin0 2 cos(* + $/2 4 
0*) x 272|yl + oP_ + |sm 2 0i(A /2 + cos0 2 ) sin0 2 cos(3* - $/2 - +o + (/)_+) 

+ 2V2|A _ | lA.+ KA^ - cos 0i ) sin 0i sin 2 2 cos(3* + $/2 - 0o_ + 4>-+) 

- 2v / 2|A 0+ ||yl + _|(A /l + cos 0i ) sin 0i sin 2 2 cos(3# + $/2 + (f> 0+ - 

- 2v / 2|A_ ||A + _| sin 2 0i(% - cos0 2 ) sin0 2 cos(3* - $/2 + <£_o - </>+-) 



A>+o ■ 
6 -- 

6-o- 



-) 



(Al) 



where TVj is the normalization constant which does not affect the angular and mass distributions. Because decays 
of vector bosons Vi — ¥ are involved, the angular distributions depend on the parameter Af i characterizing their 
decay, defined as Af = 2g v g^/(g v 2 + g^) [31]. This parameter is 1 for W decays and approximately 0.15 for 
Z — » £~£ + . Equation (|Aip represents a more general version of Eq. (Bl) from Ref. [20], where sign conventions are 
different between the two equations. Conventions for Eq. (IA1|) are consistent with Eqs. ©-([4]). The functions F/j(9*) 
are defined through the Wigner (^-functions as 7 



m=0,±l,±2 



J ( 



(A2) 



where f m are fractions of the X particle polarization as defined in Ref. [20] . In qq annihilation the resonance X can 
only be produced by m = ±1, whereas in gluon fusion m = ±2 or 0. The relative fractions of m — ±2 and are 
determined by amplitudes in Eq. (|2"Tj) which simplify in the case of couplings to two massless gluons and depend on 
production couplings in Eq. (jT5J) . The relative fraction of qq — > X production is denoted by f qq and is determined by 
the ratio of cross-sections, including effects of parton structure functions. This leads to 



/+ 



f-l = 



fzl 

2 

fz2 
~2 



fqq 
2 



■f+2 - f-2 - — = (1 - fqq)j^ 



|Af_| 2 



= 1 \ A 9 ap\ 



= (i - /,,-) 



l^+l 2 



fo - fzO - (1 - fqq) ' 



a,/3=±l 

2 + \A"_\ 2 



E 



a.3=± 



(A3) 



For a spin-zero resonance f qq = and /o = 1. For a spin-one resonance f qq = 1. For a spin- two resonance, generally 



7 The convention presented here differs from that in Ref. [20|. All probability distributions are invariant under the simultaneous 
transformations 8* — > (it — 0") and <t>i — > (-7T + "^i). The different convention is equivalent to either of these two transformations. 
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all polarizations are possible. The minimal couplings of a spin-two resonance correspond to /o = 0. Specific examples 
of Ff tj for J = 0, 1, 2 are given in Ref. 

Appendix B: Angular and mass distributions 



We illustrate MC simulation and compare it to the derived analytical angular and mass distributions in Figs. 1111 
and [T2l for the ZZ, and in Fig. [13] for the WW final states. The X — > 77 distributions are shown in Fig. [2j We have 
also validated that results presented in this paper using Eqs. (|2"3")l and (pM|) are nearly identical if in place of analytical 
parameterization of the probabilities, V, we use matrix element calculations from the vector algebra employed in the 
event generator. The two methods are conceptually independent but are mathematically equivalent, apart from the 
normalization of the probabilities which is easier to calculate with the analytical parameterization. We provide the 
necessary code for both methods in Ref. 35 [. 




m 2 [GeV] 



13 26 39 
m 2 [GeV] 



m 2 [GeV] 



m 2 [GeV] 



FIG. 11: Distributions of the observables in the X — > ZZ analysis, from left to right: spin-zero, spin-one, and spin-two signal, 
and qq — >• ZZ background. The signal hypotheses shown are (red circles), (green squares), (blue diamonds), as 
defined in Table [U Background is shown with the requirements 7712 > 10 GeV and 120 < mu < 130 GeV. The observables 
shown from top to bottom: mi and 7712 (where mi > 7712). Points show simulated events and lines show projections of analytical 
distributions. 
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FIG. 12: Distributions of the observables in the X — > ZZ analysis, from left to right: spin-zero, spin-one, and spin-two signal, 
and qq — >• ZZ background. The signal hypotheses shown are (red circles), (green squares), (blue diamonds), as 
defined in Table[U Background is shown with the requirements mi > 10 GeV and 110 < mu < 140 GeV. The observables shown 
from top to bottom: cos#*, $i, cosfi, cos 82, and $. Points show simulated events and lines show projections of analytical 
distributions. 
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FIG. 13: Distributions of the observables in the X — > WW analysis, from left to right: spin-zero, spin-one, and spin-two 
signal. The signal hypotheses shown are (red circles), (green squares), (blue diamonds), as defined in Table U 
The observables shown from top to bottom: mi,2, cos 8*, $i, cos #1,2, and Points show simulated events and lines show 
projections of analytical distributions. 
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